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Effects of spin polarization on the structure of symmetric nuclear matter and 
pure neutron matter are investigated. We show that the spin polarization 
induces a deformation of the Fermi spheres for nucleons with spin parallel 
and opposite to the polarization axes. This feature can be related to the 
CN ■ structure of the one pion exchange contribution to a realistic nucleon-nucleon 

■ interaction. While the anisotropies in the momentum distribution lower the 

. energy of the system by small amount, the associated variations of the single 

particle energies with the angle between the polarization axis and the particle 
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momentum are significant. 



The magnetic susceptibility of nuclear matter has received a considerable amount of attention 
as this quantity is relevant for the studies in astrophysics of compact objects. The strong 
magnetic fields in the interiors of neutron stars may lead to a magnetization of the material 
of nuclear matter, which implies a spin polarization of the interacting nucleons. Various 
ff~j ■ groups have discussed the possibility of a phase transition of normal neutron matter to a 
ferromagnetic state [0-01- With realistic models of the nucleon-nucleon (NN) interaction, 
however, such a phase transition seems to be suppressed for neutron matter up to densities 
well above the saturation density p of symmetric nuclear matter. Moreover, it turns out 
that the NN interaction even tends to reduce the susceptibility below the value of a non- 
interacting Fermi gas 
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where fi and m represent the magnetic moment and the mass of the neutron while fop 
stands for the Fermi momentum of system. This reduction is related to the density of states 
at the Fermi surface, which can be expressed in terms of an effective mass m*, which is 
typically smaller than the bare nucleon mass m, and the strength of the spin-spin interaction, 
represented by a Landau parameter Gq. An estimate for the susceptibility is then 

y m* , , 
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If, however, the spin-spin interaction becomes attractive at densities above the nuclear sat- 
uration density, as it is the case for some effective interactions like the Skyrme force, the 
system may undergo a phase transition to a ferromagnetic state [RBI . 
The polarizability of nuclear matter or its spin-spin correlation function is also an important 
ingredient for the study of the interaction of neutrinos with nuclear matter. This interaction 
of neutrinos with dense matter is of high importance for the simulation of supernova explo- 
sion and the cooling of neutron stars. In order to compute the mean free path of neutrinos 
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in nuclear matter one has to determine the nuclear response function for density and spin 
excitations [§-12| in particular for lower energies. These response functions are essentially 



characterized by the single-particle spectrum of the nucleons with momenta p around the 
Fermi momentum, which is often parameterized in terms of an effective mass and momentum 
independent single-particle potential U 

e<P) = ^ + U, (3) 

and the residual interaction, which is conveniently characterized by the Landau parameters 
F and Go for the density and spin response functions, respectively. 

The studies of the nuclear polarizability and the nuclear response functions all consider a 
conventional model for nuclear matter, where the single-particle states are occupied accord- 
ing to the magnitude of the single-particle momentum p. As a typical example we refer 
to the Brueckner-Hartree-Fock (BHF) calculation of Vidana et al. f7|, who evaluate the 
single-particle energies for nucleons with momentum p and spin projection s according to 

p^ r (Pq 
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where G-matrix is the solution of the Bethe-Goldstone equation for a realistic model of the 
NN interaction, u is the starting energy and n(q, s') is the occupation probability of nucleons 
with momentum q and spin projection s'. In the BHF approximation for systems at tem- 
perature T = this occupation probability is represented by a step function which restricts 
the occupation to states with momenta below the Fermi momentum kp(s') of nucleons with 
spin projection s'. Note that in order to shorten the notation, we have suppressed all labels 
referring to the isospin of the nucleons. 

It is well known that a realistic NN interaction contains strong tensor and spin-orbit com- 
ponents, which couple the spin and momentum or spatial dependence of the interaction. In 
the case of the deuteron these tensor components lead to dumbbell shape if one considers 
deuterons with spin projection M s = ±1 flT3|. These considerations lead to the conjecture 



that also the momentum distribution of polarized nuclear matter could be deformed with 
respect to the axis of spin polarization. Configurations with deformed momentum distribu- 
tions could be lower in energy than the spherical ones discussed so far, which would affect 
the polarizability. A deformed momentum distribution would also be connected with the 
single-particle spectrum, in which the single-particle energies depend on the magnitude and 
the direction of the momentum p of the nucleon under consideration. This would lead to 
different response functions and influence the mean free path of neutrinos as well as other 
properties of the nuclear matter. 

In order to investigate this question, we consider a Fermi momentum that depends on the 
angle $ with respect to the spin-polarization axis, parameterized in terms of 

k F (s, <&) = k F0 (s) (l + 5(a) cos 2 tf) , (5) 

a parameterization which turned out to be useful. With this assumption we can calculate 
the contribution of the single-particle potential U(p,s) to the single-particle energy in (J|), 
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employing the partial wave expansion for the G-matrix as outlined e.g. by Haftel and Tabakin 




U(p,s)= J2 I jf^n(q, s') (L'M L ,Ss + s'\JM)(LM L ,Ss + s'\JM) 
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with K cm and k r referring to the center of mass and relative momentum of the interacting 
pair of nucleons, respectively. The brackets (LrriL,S, M$\JM) denote the Clebsch-Gordan 
coefficient for the coupling of the spin and orbital angular momentum of the relative motion. 
Note that for a system without spin-polarization the spherical harmonics Yl,m drop out. If, 
however, we assume a spin-polarization Eq. (||) yields a single-particle potential which 
depends not only on the magnitude of the momentum p but also on the angle relative to 
the axis of spin-polarization, even if we assume spherical momentum distributions n(q,s'), 
i.e. S(s') = in fl5|). From the momentum dependence of the resulting single-particle energies, 
one can find a Fermi momentum according to (|5|) and recalculate the single-particle energy 
with this deformed momentum distribution until a self-consistent solution is obtained. 
Results for the single-particle potential of such BHF calculations are displayed in Fig. || for 
nuclear matter at the empirical saturation density po with a spin polarization 

Ap = ^ (7) 
Po 

of 25 percent. The calculation have been performed using the charge-dependent CDBonn 



NN interaction of Machleidt et al. ||15|| . The spin asymmetry as well as the angle-dependence 
of the single-particle spectrum have been ignored in the Pauli operator and starting energy, 
which were used to solve the Bethe-Goldstone equation. The single-particle potential for the 
nucleons with spin parallel to the polarization axis (left panel of Fig. |1[) and a momentum 
of 1.46 fm _1 (the Fermi momentum of the spin up momentum distribution if it would be 
spherical) is about 2.7 MeV more attractive in the direction of the momentum perpendicular 
to the polarization axis (6 = 7r/2) than for a momentum parallel or antiparallel to this axis. 
This means that we obtain an oblate momentum distribution for the nucleons with spin 
parallel to the polarization axis, which can be parameterized according to Eq. (|5]) with 
a deformation parameter 5(f) « —0.03. It is worth noting that this oblate deformation 
is in line with the observation made for the deuteron |13| as the dumbbell shape of the 
density distribution mentioned above corresponds to a momentum distribution with oblate 
deformation. 

The situation is opposite for the nucleons with spin antiparallel to the polarization axis. In 
this case we observe a smaller single-particle energy for momenta parallel to the polarization 
axis as compared to the momentum perpendicular to this axis (see right panel of Fig. p. 
This leads to a prolate deformation of the momentum distribution of the nucleon with spin 
antiparallel to the polarization axis. 

The angle-dependence of the single-particle potential is of course much stronger in the case 
of totally polarized nuclear matter. In this case we obtain a energy difference of around 
20 MeV between the single-particle energies for the direction of the momentum parallel as 
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compared to perpendicular to the polarization axis, if we consider a momentum of 1.72 
fm _1 , which is the Fermi momentum of the spherical distribution for this polarization. This 
corresponds to a deformation parameter (5(f) m —0.15 in the parameterization of Eq. (H). 
The effects are weaker for pure neutron matter. For example, in the case of completely 
polarized neutron matter at density po, the gain in the single-particle energy is around 
6 MeV if the momentum is perpendicular to the polarization axis z as compared to the 
same momentum parallel to z (assuming a neutron with spin parallel to z and a momentum 
p — 2.16 fm _1 which is the Fermi momentum of the spherical distribution). This corresponds 
to a deformation parameter 5(f) « —0.03. 

To understand the origin of this deformation we recall that within the Hartree-Fock ap- 
proximation the non-locality or momentum dependence of the self-energy of a nucleon in 
nuclear matter is mainly due to the Fock-exchange contribution in (J3j). For the sake of 
simplicity let us consider completely polarized nuclear matter. In this case the contribution 
of the one-pion exchange to the Fock term in Eq. (|]) for a nucleon with spin parallel to the 
polarization axis can be described by 



[&x ■ k)(a 2 ■ k) 



m 2 + k 2 



gf,pf )n(g,f). (8) 



In this equation we have approximated the one-pion exchange term by its static, nonrela- 
tivistic reduction in which ex, represents the vector of Pauli matrices acting on nucleon i and 
k denotes the momentum transfer which is given in terms of the momenta of the interacting 
nucleons by k = p — q. Since the spin part of the matrix element in (|8|) yields results different 
from zero only for the ^-component of the spin matrices (z-axis parallel to the polarization 
axis) one obtains 

At4(£f) ~ / d 3 q iPz ~ q ^ n(q,]). (9) 



m 2 + k 2 

Indeed, we find that the self-energy for completely spin-polarized nuclear matter increases 
with the absolute value of p z , but is rather weakly dependent on p x and p y . 
The situation is different, if we consider the self-energy of a nucleon with spin antiparallel 
to the axis of spin polarization. In this case we have to consider matrix elements of the form 



d 3 q (P f, q f 
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k 2 
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which lead to a self-energy which increases with the component of p perpendicular to the 
z-axis. This means that the angle-dependence of the self-energy can be understood qualita- 
tively in terms of the Fock contribution of the pion-exchange part of the NN interaction. 
We will turn next to the effect of the deformed momentum distribution on the total energy. 
The energy per nucleon is given by 

E = E s Jd 3 pn(p,s)±(^+e(p,s)) 

A~ Z s fd 3 pn(p,s) ■ [ > 
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One finds that the gain in energy due to the deformation of the momentum distribution 
is rather small (see table |). Even in the case of completely polarized nuclear matter, the 
deformed momentum distribution yields an energy per nucleon which is 0.53 MeV smaller 
than the energy calculated for the spherical momentum distribution. This is a rather small 
correction compared to the energy difference between the polarized and unpolarized energy. 
Therefore we cannot expect that the effects of deformed momentum distributions will alter 
the conclusions of e.g. J7J about the equation of state for spin polarized matter or the 
possible existence of a ferromagnetic state of nuclear or neutron matter. The situation, 
however, might be compared to the deformation degrees of freedom in finite nuclei. The 
possibility to form deformed nuclei is not very relevant for the total energy of finite nuclei. 
The gain in energy due to deformation of nuclei is typically of the order of 0.05 MeV per 
nucleon or even less (see e.g. [|HJ). The deformation, however, has rather significant effects 
on the excitation spectrum at low energies, which is quite different for spherical as compared 
to deformed nuclei. Also, the dependence of the single-particle energies on the direction of 
the nucleon momentum in spin polarized nuclear matter may lead to quite different response 
functions at low excitation energies. This may have a significant effect e.g. on the neutrino 
propagation in neutron stars with large magnetic fields. 
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unpolarized 




completely polarized 








spherical 


non-sph. 


*(T) 


Nuclear Matter 


-16.46 MeV 


5.71 MeV 


5.17 MeV 


-0.12 


Neutron Matter 


11.58 MeV 


56.94 MeV 


56.91 MeV 


-0.03 



TABLE I. Energy per nucleon for nuclear matter and neutron matter of density po = 0.17 fm . 
Results are presented for the unpolarized system and the system with complete spin polarization. 
For this case we show the energies assuming a spherical and non-spherical momentum distribution. 
The last column shows the deformation parameter for the self-consistent solution according to Eq. 
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FIG. 1. Single-particle potentials for nucleons with s pin parallel (left panel) and antiparallel 
(right panel) to the spin polarization function of the angle 9 between their momentum 

and the polarization axis in nuclear matter. Results are given for various values for the modulus 
of the momentum p (expressed in fm -1 as indicated in the legend). The single-particle potentials 
are symmetric around 9 = n/2. 
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